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Abstract. 

We construct a closed-form, fully analytical 4-metric that approximately represents 
the spacetime evolution of non-precessing, spinning black hole binaries from infinite 
separations up to a few orbits prior to merger. We employ the technique of asymptotic 
matching to join a perturbed Kerr metric in the neighborhood of each spinning black 
hole to a near-zone, post-Newtonian metric farther out. The latter is already naturally 
matched to a far-zone, post-Minkowskian metric that accounts for full temporal 
retardation. The result is a 4-metric that is approximately valid everywhere in space 
and in a small bundle of spatial hypersurfaces. We here restrict our attention to quasi- 
circular orbits, but the method is valid for any orbital motion or physical scenario, 
provided an overlapping region of validity or buffer zone exists. A simple extension 
of such a metric will allow for future studies of the accretion disk and jet dynamics 
around spinning back hole binaries. 
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1. Introduction 

A closed-form, analytic understanding of the two-body problem in General Relativity 
has proven quite elusive. From the point of view of gravity, this holds true for Newtonian 
stars described by Newtonian gravity, but even more so for binary compact objects, 
such as black holes (BHs) and neutron stars (NSs), in the last stages of inspiral and 
merger. Perhaps, this can be traced to the fact that the "two-body problem" in General 
Relativity is really a "one-spacetime" problem, where the orbital dynamics cannot 
be easily separated from the spacetime dynamics. Fortunately, the breakthroughs in 
numerical relativity [U [2J E] have now made possible to fully compute the dynamics of 
any binary BH (BBH) spacetime for a wide variety of mass ratios and spins parameters. 
However, these calculations are still very expensive, and therefore are limited to few 
orbital evolutions with the BHs starting at relatively close separations. 

A closed-form, analytic understanding of the full spacetime in the inspiral regime 
would be extremely useful to study a variety of astrophysical phenomena that require a 
large number of orbits. The study of the dynamics of accretion disks around BBHs and 
of any associated electromagnetic emission and jets is a good example. With this goal 
in mind, we develop an asymptotically matched spacetime metric for non-precessing, 
spinning BBHs. Further motivation for our work is described in [3]. 

Our work builds on important previous analytic perturbative techniques, which 
have been successfully developed to tackle the BBH problem. Post-Newtonian (PN) 
theory (see, e.g., [5] for a review) can be successfully used to describe the motion of 
the BBH in the early stages of the inspiral. Here, all fields are perturbatively expanded 
in a slow motion !)/c < 1 and weak field GM/(rc 2 ) <C 1 approximation [|J Close to 
each of the BHs, one can use perturbation theory, where all fields are treated as small 
deformations of a known analytic solution, such as the Schwarzschild or Kerr metric. 

During the late stages of a binary BH, however, the whole spacetime cannot be 
described using a single perturbative method. For clarity, let us classify different spatial 
regions on a spacelike hypersurface into zones [B] (see Figure [I] and Table [I] for a 
BBH). The inner zone (IZ) is defined as the region close enough to either BH that 
the metric can be treated as a perturbation of the Kerr spacetime due to some external 
universe El El [TO], HH [12]. The near zone (NZ) is defined as the region both far 
enough from either BH that the metric can be perturbatively expanded in a weak-field 
approximation, while simultaneously much smaller than a gravitational wavelength, so 
that time retardation can be treated perturbatively [5]. The far zone (FZ) is defined 
as the region sufficiently far away from the center of mass of the system (much farther 
than a gravitational wavelength) that the metric can be modeled via a multipolar post- 
Minkowskian expansion [5J. In each of these zones, one can construct an approximate 

| Here, c and G are the speed of light and Newton's gravitational constant, while v, M and r are 
the characteristic velocity, mass and size or separation of the system. Henceforth, we will adopt the 
geometric unit system, where G = c = 1, with the useful conversion factor 1M Q = 1.477 km = 
4.926 x 1CT 6 s 
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4-metric, in some coordinate system well-adapted to that particular zone, and in terms 
of certain physical parameters, like the BH masses and spins. 




Figure 1. BH1 and BH2 are denoted by solid black dots, where the orbital separation 
is b. The BZs are denoted with cyan shells, the outer one representing the FZ/NZ BZ 
and the two inner ones representing the NZ/IZ BZs (see also Table [T]). The circular 
nature of these buffers zones is because the diagram is schematic; in practice, one 
expects these to be distorted. The IZ, NZ and FZ are also shown in the figure. 

Asymptotic matching is the mathematical technique that allows one to relate the 
coordinates and parameters used in adjacent zones inside of a common overlapping 
region of validity, or buffer zone (BZ) (see Figure [l] and Table [TJ. Formally, asymptotic 
matching requires that one set the asymptotic expansion of approximate 4-metrics inside 
the BZ equal to each other. This yields a system of differential and algebraic equations 
for the coordinate and parameter transformation that relates adjacent approximate 
metrics. After applying such a transformation, one can then join the approximate 
metrics with carefully constructed transition functions to obtain an approximate global 
metric. 

Such an approximate global metric is ideal as initial data, because in practice 
asymptotic matching is carried in the neighborhood of some fiducial time to, i.e., the BZ 
is the product of a 3-torus with a small segment of the real (temporal) line. Alvi [T31 IT4] 
was the first to attempt such an initial data construction, but ended up carrying out 
asymptotic patching rather than matching [§] Yunes et al. [151 EEO US HH] succeeded in 
carrying out matching for non-spinning BBHs and this data has recently been evolved 
in [20] (see also [21] for numerical evolutions of superposed tidally perturbed BHs) . 

In this paper, we extend the calculation in [TS1 [TTJ, [12] to non-precessing, spinning 

§ When patching, one sets the metrics equal to each other at a point, instead of in an entire BZ region 
(for more details see [131 Ho] V 
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Table 1. r; n and r out denote the approximate inner and outer boundary radius, 
respectively, while "Region" stands for the spatial domain in which the BZs exist, ta 
is the distance from a ^4th BH with mass m A , r is the distance from the center of mass, 
and b and A are the orbital separation and the gravitational wavelength, respectively. 



Zone 


^in 


r out 


Region 


IZ BH1 (n) 









IZ BH2 (r 2 ) 





< b 




NZ (r A ) 


> m A 


< A 




FZ (r) 


>6 


oo 




IZ-NZ BZ 






771^4 <C r A b 


NZ-FZ BZ 






b < r < A 



BBHs in a quasi-circular orbit, with BH spins aligned or anti- aligned with the 
orbital angular momentum. The NZ and FZ are modeled via the standard PN 
approximation [22j |23| I24"t |25| [26] (and references therein). The IZ around each BH 
is modeled as a Kerr spacetime with a leading-order vacuum perturbation, following 
Yunes and Gonzalez [27] . The IZ metrics are cast in a coordinate system that is both 
horizon-penetrating and harmonic [2E]- This simplifies the matching and allows for 
excision of the singularities. 

We first carry out asymptotic matching in the non-spinning limit. Our results for 
the matching coordinate and parameter transformation differ from those of pj)J [TTJ, H~9] 
because our IZ metrics differ by a gauge transformation. We then carry out asymptotic 
matching for an aligned or counter-aligned, arbitrarily spinning binary. The matching 
coordinate and parameter transformation are not affected by spin to the matching order 
studied. Asymptotically matching the metrics to higher order would introduce spin- 
dependent terms in the coordinate and parameter transformation, but this would require 
perturbed Kerr metrics valid to higher multipolar order in the IZs. Once asymptotic 
matching has been carried out, we stitch the different metrics via appropriate transition 
functions, thus obtaining a global approximate metric. 

We verified that this metric is indeed an approximate solution to the Einstein 
equations. One might be worried that the stitching procedure introduced errors 
in the global metric that are larger than those inherently contained in any of the 
approximate metrics. This is not the case because the transition functions used satisfy 
the Frankenstein theorems of [18]. We verified this qualitatively by visually inspecting 
different components and the volume element of the global metric. We then verified 
this quantitatively by evaluating the Ricci scalar for the global metric on a t = 
spatial hypersurface. We find that the global metric is an approximate solution to 
the Einstein equations for all values of spin considered, provided the binary orbital 
separation is sufficiently large, so that BZs in which to carry out asymptotic matching 
exist. We further verified that as this orbital separation is increased, the satisfaction of 
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the Einstein equation increases at the expected rate (given by the matching order). 

Such a global metric is technically valid close to a t — spatial hypersurface 
because the BZ in which matching is carried out is formally a small 4-volume about 
this hypersurface. One can, however, extend this global metric to a dynamical global 
representation of the full spacetime, following the procedure in [29J. This scheme 
essentially consists of using a sequence of temporally-spaced global metrics and properly 
gluing these together. Such a spacetime was successfully employed to model how the 
inspiral of BBHs affects accretion disks [I] . 

The remainder of this paper is organized as follows. In section [2j we summarize how 
the metrics of each zone are constructed. In section [3j we first redo asymptotic matching 
between IZ and NZ metrics in the non-spinning case but using an IZ metric different 
from the one used in [TSJ [TTJ, [TBI US] > an d then extend it to the spinning case. We present 



some technical details and supplementary analyses in the appendices. In Appendix A 



we review how this paper is related to |27j, and extend the formulation by including the 
m = azimuthal mode. In Appendix B we present a detailed and complete analysis of 
asymptotic matching. In Appendix C , we discuss how to glue the metrics in each zone 
with transition functions. In Appendix D, we summarize some higher order metrics 
which are available for the non-spinning case in the supplemental material of [T9] . 

All throughout, we use the following conventions, following mostly Misner, Thorne 
and Wheeler [30]. We use the Greek letters (a, (3, • • •) to denote spacetime indices, and 
Latin letters • • •) to denote spatial indices. The metric is denoted g^ u and it has 
signature (— , +, +, +). We use geometric units, with G = c = 1. 



2. Approximate Metrics 

2.1. Inner Zone 

The metric in either IZ is approximated with a Kerr solution plus a linear vacuum 
perturbation [2"7j : 

C = 917 + fcj ; (i) 

the construction of h% v is difficult but essential for our purposes. Vacuum perturbations 
of Schwarzschild have been studied, for example, in |31j, where one can use the 
Regge-Wheeler-Zerilli-Moncrief formalism [32| 133} IM] . When considering the Kerr 
background, however, the perturbation equations do not easily decouple, as tensor 
spherical harmonics are not eigenfunctions of the angular sector of the linearized Einstein 
equations [35]. Instead one has to carry out the so-called Chrzanowski procedure |36j, 
amended by Wald [37] and by Kegeles and Cohen [38], as we describe below. 

In the Chrzanowski procedure jSU], the vacuum perturbation is obtained from 
a so-called Hertz potential \I/ via 

K v = V [*] , (2) 

where [•] is a differential operator. This potential must satisfy a certain non-linear 
differential equation with a source given by the Newman-Penrose scalar ip (or t/> 4 
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depending on the radiation gauge used). This differential relation can be inverted 
through the Teukolsky-Starobinski relation to yield ^ in terms of ^0,4 ([311 EH HO]). 
Thus, the construction of a metric perturbation reduces to finding an appropriate 
solution for the Newman-Penrose scalar. 

The Newman-Penrose scalar must, of course, satisfy the Teukolsky equation, 
but when the source is a slowly-varying external universe (as is the case when the 
perturbation is due to a BH binary companion in a circular orbit with large semi- 
major axis), this equation can be solved perturbatively. One first performs a harmonic 
decomposition of ipo in terms of spin-weight 2 spherical harmonics (or ip^ in terms of 
spin- weight —2 spherical harmonics), ^Y 1 ™". 

^o = J2 R im(r)z im (v) 2 Y £m (6, 0) , (3) 

l.m 

where the radial and temporal dependence are product decomposed in terms of unknown 
real functions Ri m (r) and unknown complex functions Z£ m (v), where v is the advanced 
Kerr-Schild time coordinate. The latter can be written in terms of certain combinations 
of electric and magnetic tidal tensor components, which characterize the perturbations of 
the external universe. To leading order in a small-hole/slow-motion approximation [JT] 
in BH perturbation theory, it suffices to keep only the £ = 2 quadrupolar deformation 
(see also Appendix A). The radial functions are then required to satisfy the (time- 
independent) Teukolsky equation, which one can solve in terms of hypergeometric 
functions [HJ [27]. This means that zg m is constant or slowly varying in time. With 
the Newman-Penrose scalars under control, one then proceeds to compute the Hertz 
potential and from that, the metric perturbation. Reference [27] provides explicit 
expressions for g™ in Boyer-Lindquist (BL) coordinates. 

For ease during the matching procedure, one wishes to express the IZ metric in a 
coordinate system that closely parallels that used in the NZ. We thus transform the 
IZ metric from BL coordinates (£bl, r BL, #bl, 0bl) to harmonic, horizon-penetrating 
coordinates (T, X, Y, Z). Harmonic coordinates are those that satisfy Dx M = 0, which 
of course defines a class of coordinate systems. Several coordinate transformations 
between BL and different members of this class were developed in [121 [53], [28], HU 05] H6J 
HTj : we here employ one such transformation that leads to harmonic coordinates that 
are also horizon-penetrating [28], namely 
rl + A 2 , r BL - r + „ 

V BL - M) COS 6>bl 



T = t BL + -± In 

r _ r _ 



r BL - r. 



(4) 



X + iY = (r BL -M + i A)e**« sin 6 BL ; IK = 0bl + - In 

r + - r_ r B L 

where A = S/M is the dimensional Kerr spin parameter associated with the Kerr 
background with spin S, and M is the mass of the Kerr background, with r± = 
M ± a/ M 2 — A 2 the location of the unperturbed inner and outer horizons in BL 
coordinates. The quantity 4>ik should be thought of as an azimuthal variable of ingoing 
Kerr (IK) coordinates, not as a spherical polar $ coordinate associated with harmonic, 
horizon-penetrating coordinates. 
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The IZ metric is then in harmonic, horizon-penetrating coordinates X a and it 
is characterized by the parameters A a = (M, A, z^ m , zi tm ), where the first two are 
associated with the background and the latter two are related to the metric perturbation, 
i.e., the real and imaginary parts of z 2m . In section |3j we will carry out asymptotic 
matching and relate these coordinates and parameters to the NZ ones. 



2.2. Near Zone 

The NZ metric is chosen to be given by the PN expansion, 

C = + Ki , (5) 

where rj^ is the Minkowski metric, while is a PN metric perturbation. The latter can 
be decomposed into spin-independent and spin-dependent terms. All spin-independent 
terms are given explicitly in [23] , and we will here only keep the first non- vanishing spin 
terms in the metric perturbation, which can be found in [2H [25] . 

The metric perturbation of linear-in-spin contributions can be written as 

6hZ' {S) = 25V^ + 0(v 7 ) , 5C (S) = -m (S) + 0(v 6 ) , 6hl z ' {s) = 0(v 5 ) . (6) 

The PN ordering system we employ is as follows: a term that corrects the leading-order 
coefficient of some expression by a quantity of relative 0(v 2n ) = 0(v 2n /c 2n ) is said to 
be of ra-th PN order. The spin potentials to leading order are 

6V<® = \e ijk v\s{n\ + (1 2) + 0(v 5 ) , 

5VP = ^e ijk s{n\ + (1 2) + 0(v 3 ) . (7) 

Here, s\ denotes the spin angular momentum of the Ath PN particle, which has 
dimensions of (mass) 2 , while the Ath particle's mass, location and velocity are given 
by niA, y\ and v\ respectively. We also use the notation, = |x — y^l and 

n A = - V l A)l r A- 

Putting all of this together, one has the 1.5PN order NZ metric, 



m 7 -, 2m 1 rti\ r . , , , 91 m 2 
flSo + 1 = -zr + 7 1 i^i ~ ( n i • vi) 2 ] - 2-i ~ m x m 2 
' i ' i 'i 



2 n r? 



rir 2 2b 3 2r 2 b 3 2r 2 b 
+ ^^(n 12 • v 12 ) + ^vlsinl + (1H2) + 0(v G ) , 



gg = - - ^e ijk 4n1 + (1 o 2) + 0(v 5 ) 

7*1 ri 



2m 



NZ - = — ^ + (1 2) + 0(v') , (8) 



ri 

where we have introduced the notation 



6 = |yi — Yal , ni2 = (yi - Y2)/b, v 12 = vi-v 2 . (9) 

The quantity b is the same as the commonly used ri 2 in the PN literature. Quadratic 
spin term are here neglected, as they enter at higher PN order. The spin-independent 
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terms are available up to 2.5PN order from [23], and implemented as supplemental 



material in [19]. We briefly summarize these higher order terms in Appendix D 



For quasi-circular orbits, certain simplifications are possible. First, we have 
ni2 ■ Vi j 2,i2 = 0(v 5 ), since the orbit radially decays due to gravitational radiation 
reaction. Ignoring the latter, b = b {cos ut, sinotf, 0}, where 



to 



^(^ <«» 

is the orbital angular frequency. Here, m = mi + m 2 is the total mass of the system. In 
section [3j we will use the positions of the holes, 



yi 



m 2 
m 



b + O^), y 2 



m 



:ir 



that are obtained in the calculation of the center of mass of the system. Higher-order 



expressions for these quantities are given later in Eq. (15). 



The NZ metric is then given in harmonic PN coordinates x a , and characterized 
by the parameters X a = (mi, 777.2, b, s\, s\). In section [3J we will carry out asymptotic 
matching and relate these coordinates and parameters to the IZ ones. 



2.3. Far Zone 



The FZ metric can be obtained via the direct integration of the relaxed Einstein 
equations H9"j I5D] in harmonic gauge, or alternatively via the multipolar formalism 
of Blanchet, Damour, and Iyer (BDI) [EO E2 E31 E]. The FZ metric in @8j [50j CL9] is 
given by 

1 



9ll 



9ok 



i-lKl + l K) 2 



-hf z + 0(v" 



9lt 



l + -h 00 



5 kl + 0(v 4 



(12) 



where the metric potentials h^z are 



h 00 

"fz 



4- + 29 H 
r 



rkl, 



u 



- 2 -B 

g u klm 



T 



Mm 1 



U 



+ 7 ^ + ^ 



h 



kl 



-2d l 
X kl (u 



Z kt (u) 



3 p 



l klp (u) 
r 



_° e ps(k\ d 
3 



j p V(u) 
r 



I 2 



+ ^h k fi + 0(v ti 



(13) 



with r = |x| the distance from the binary's center-of-mass to the field point, n k := x k /r, 
and u = t — r is the retarded time. We follow here the PN order counting explained 
in [19], and we keep terms up to the same order as those retained in the NZ metric. 
Higher-order terms are of course available, and implemented for the non-spinning case 



as the supplemental material in [19] (see Appendix D for a brief summary). 
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The FZ metric is then given in terms of source multipole moments, which are 
corrected by spin terms. The spin contributions to the multipole moments can be read 
from Eqs. (B5) and (CI) in [M] (see also [26]): 

1 



6J, k 



(S) 



S l S 2 i 



5J ( kl 



(S) 



(3 s\y[ - 6 kl s?y™) + (1 2) 



j a and s l A are the same as in the NZ but evaluated at retarded time u. 



(14) 
All 



where y\, v 
non-spinning terms are given explicitly in [19J. 

We have verified that the NZ and FZ metrics are automatically asymptotically 
matched in the BZ constructed from the intersection of the NZ and FZ. In practice, 
the leading 1.5PN order spin effects in the NZ metric, g^o an d 9oi are matched to the 
spinning parts of (1/2)/^ and — h® k in the FZ metric and g$f, respectively. That is, if 
one asymptotically expands both of these metrics in this BZ, then they are automatically 
equal to each other without requiring any coordinate or parameter transformation. 

To establish this result, one need to only worry about the following transformation 
of the relative to center-of-mass coordinates. According to Eq. (5.5) of [25], the relative 
position b n 12 = yi — y 2 and velocity v 12 = v x — v 2 are converted from the quantities 
in the center of mass system via 



yi 



y2 



m 2 7i 

m 2 
mi ri . 

— + i$ 

777 2 



V 



v 2 - 



777 V 

j). 

777 

b 



bn 



V 



12 



711 



V12 X £ 



bn 



7/ 



12 



III 



Vi 2 X S , 



(15) 



up to 1.5PN order. Here, 6 = (mi— 777 2 )/777, rj = 777 i777 2 /m 2 and S = 777(S 2 /777 2 — Si/mi). 

The above transformation is affected by the choice of spin supplementary condition 
(SSC) (see, e.g., Appendix A of [551 El])- Although the NZ and FZ metrics are 
technically computed using different SSCs, this is not a problem here due to the PN 
order to which we work. 



3. Asymptotic Matching 

In this section, we carry out asymptotic matching between IZ and NZ metrics. We will 
here follow the same procedure as that first introduced in [15], and recently refined 
in [TH]. We first concentrate on matching the IZ metric around BH1 and the NZ 
metric; matching between the other IZ metric and the NZ can be obtained later via 
a symmetry transformation. In asymptotic matching, one expands both metrics in 
the BZ, m\ < ri < 6 and t < 6, and then requires that they be diffeomorphic to 
each other. This leads to a set of differential equations that relate the coordinates 
used in each metric, as well as a set of algebraic equations that relate the parameters 
used in each zone. The IZ metric around BH1 depends on 5 complex parameters z 2m 
(m = (—2, —1, 0, 1, 2)), which must be determined with the matching procedure. 

Before proceeding with the matching, let us briefly summarize the coordinates, 
parameters and expansions which we employ on the IZ and NZ metrics. The NZ 
metric is expressed in NZ, harmonic, PN coordinates x a and depends on parameters 
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A" = (mi, m2, b, s\, s 2 ), where we recall that mi,2 are the masses of the PN particles, 
s\ 2 is their spin angular momentum, and b is their separation. The IZ metric is expressed 
in IZ, harmonic, horizon penetrating coordinates X a and depends on parameters 
A a = (M, S, ZR !m , zi <m ), where M and S are the mass and the magnitude of the spin 
angular momentum associated with the Kerr background, while zj^ m and zj m are the 
real and imaginary parts of the 10 tidal field parameters that characterize the metric 
perturbation. We expand these metrics in the BZ in powers of (m 2 /b) 1 ^ 2 = 0(y): 



' /9 



8=0 
71 



' /2 

A ° = E (x) m * OT + ' (16) 

or more explicitly 

n 72 

M ( Aa )=E(?y ( M ).M + °K +1 )' ( 17 ) 



o V b } 

S ^ a ) = E (t) (a/3) + ' (18) 

i=0 

^//, m (A Q ) = e (x) M + • (19) 

i=0 

Here, {X a )i are z functions of the NZ coordinates x@ , while (A a ) i are i functions of the 
NZ parameters A' 3 . Similarly, (M)i, (S)i and (zR/i >m )i are also functions of the X 13 NZ 
parameters. We here take the sums up to n = 2, i.e., carry out asymptotic matching to 
0[{m 2 /bY}. 

With this expansion in hand, we proceed to carry out asymptotic matching in the 
next subsections. We begin by focusing on the non-spinning case. This is different from 
the work in [15l [19] because we here use the IZ metric of [27] in the non-spinning limit. 
We will then proceed with the matching of the spinning case. 

3.1. Expansion of the Non-Spinning IZ and NZ Metrics 

Let us first expand the NZ and IZ metrics in the BZ. Formally, we expand the NZ metric 
as in [19]: 

9«p = (<M)o + J^-{g*p)i + (x) (9ap)2 + 0(v 3 ) , (20) 
without any spin contributions, where 

Wo = (s5)i = o, 



la}) 2 



2m 1 b 2 , 

' 2 - - i r 



i)o ■ (b)o} + \ 2 {S[(ri) ■ (b) ] 2 - [(n) ] 2 }] A Q/3 . (21) 



m 2 (ri) b I J 6 2 

We have here defined a "lowered four dimensional Kronecker delta", A a p = 
diag(l, 1, 1, 1), and (6 fe ) = x k = {1, 0, 0} is a unit vector, with the assumption that each 
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BH is initially located on the x-axis. We also use the Cartesian PN coordinate basis 
vectors, i a = {1, 0, 0, 0}, x a = {0, 1, 0, 0}, y a = {0, 0, 1, 0} and z a = {0, 0, 0, 1} 
later. 

Before we expand the IZ metric in the BZ, it is convenient to re-express the 
parameters z R ^ m and z^ m in terms of the electric Ski and magnetic B^ tidal tensor 
components: 

Z R,0 — ~ 2 &XX — 2 Syy , Zjfl = —2 Bxx ~ 2 Eyy , 

= -2 - 2 £y Z , = -2 £ XZ + 2 £ yz , 

zj,i = 2 £ yz - 2 £ xz , = -2 £yz - 2 , 
^,2 = -2 (^xx - ^yy) - 4 i3 X y , zr—i = -2 - £yy) + 4 i3 X y , 
z/, 2 =4S XY + 2 {B YY -Bxx), zi,-2 = -4 £xy + 2 (Syy -Bxx), (22) 
where we have used the traceless conditions, 

Sxx + S YY + £zz = 0, B xx + B YY + B ZZ = 0, (23) 



(see also Appendix A). This will allow for a more direct comparison of our calculation 
with those of [T9] . 

Let us now expand the IZ metric in the BZ. Since the full form of the BZ expanded 
IZ metric is too long and unilluminating to be included here, we only provide terms up 
to the second order: 

(9^)o = V«p, (5^)1 = 0, 

Wolh = 2 -^^y Q - p(^)o(A fe ) (X% , 

(flSDa = ^^ {£klUXk)o{Xl) ° + W {R)oiEik)o{Xk)o ' 

= (^Wo ~ 3P ( ^ )o(X * )o(X ' )o )^' " W^ lMR)l ■ (24) 
Recalling here that the tidal tensors must also be expanded as in Eq. ( Jl~9| , we have 

£h = ^(S«)o + 0(v 3 ) , B kl = (^) 3/2 i(g„) + (9(t; 4 ) , (25) 

to leading order (see also Eq. (5.3) in [E]). Since is higher order than E^i, it 
can be ignored when matching to leading order. This metric is nonsingular at the 
horizon R = M due to the use of harmonic, horizon-penetrating coordinates; this is also 
true in the spinning case. Therefore, for sufficiently small perturbations, the horizon- 
penetrating character of the metric is preserved. 

3.2. Matching in the Non-spinning Case 

We carry out asymptotic matching order by order in (1112/ b) ly/2 and use the same notation 
as in [19] . 



Asymptotically Matched Spacetime Metric 



12 



3.2.1. Zeroth-Order Matching: 0[(m2/b) ] At zeroth order, we have 

«)o = (A* 7 )o( V)o(s5)» > ( 26 ) 
with = d a X^ . Because (g^o — (9ap)o — Va/3, the matching condition reduces 
exactly to that of [19]. Taking into account the position of BH1, we then have 

(X a ) = x a - — b x a = x a . (27) 
m 

3.2.2. First-Order: 0[(m2/b) 1 ^ 2 ] At first order, we have 

G$)i = (A» 7 )o(V)o(»5)i + 2 (Vlif^lloWlo , (28) 
and using {g™f)i = (fl^i = °> (A/)o = 8j, and (g™ 5 ) = r} a/3 , the above equation 
becomes 

(A M ))i = 0- (29) 
Equation (5.6) of [TJJ] can be used to find the general expression of (X a )i. 

(X a )i = (F a/9 )x^ + (D a )i = (F Qj8 )i^ + (C a )x , (30) 

where (-F aj 3)i is a constant 4x4 antisymetric matrix, (C a )i (also (-D a )i) is a constant 



4x1 matrix and we have used Eq. (27). This is the most general solution of the flat-space 
Killing equation. 

3.2.3. Second-Order Matching: 0[(m,2/b) 1 ] At second order, using (g^aji = 0, we have 



(g%) 2 = (A^) (A fi %(g™ 6 ) 2 + (V)i(V)iW)o + 2(A ( ^) 2 ( V)oG#)o , (31) 
which gives us 

«) 2 = {9%)2 + (F Q 7 )i(^ 7 )i + 2{A Wi )* . (32) 
When we define 

S a p = A( a (3) = d( a Xp) , (33) 



Eq. (32) can be rewritten as 

2(5^) a = (fl5) a - (<&) 2 - (^a 7 )l(^)l , (34) 
which is the matching equation one must solve at second-order. 

Based on [TH], the above equation is integrable if the Riemann tensor associated 
with S a p is null at second order, i.e., 

Za0y8 = d a p{S~ i $)2 + d 1 s{S a p)2 — d a s(S^)2 — d^(S a s)2 = , (35) 

for all sets of a, f3, 7 and 5. When we take 7 = 5 = and a, (3 = i,j 7^ 0, the 
integrability condition becomes 

<%(#oo) 2 = 0. (36) 

By linear independence, we can split this equation into a polynomial and a non- 
polynomial part. The latter (the one that diverges when — > 0) gives 
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where we have introduced the notation (rf)o = x from the zeroth-order matching. 
This equation then forces (M) = mi. For the polynomial part, it is easy to verify 
that constant and first order pieces (in the sense of a polynomial decomposition 
a + bX + cX 2 + ...) give consistent, although trivial, equations. The quadrupolar piece 
gives 

] 2 - [( ri ) ] 2 + (E kl ) (x k Ux l ) 



d. 



'j 



3[(ri)o-(b) t 







(3* 



which then forces 

\^ijjO &ij 3XiXj . 



(39) 



This is consistent with Eq. (5.17) of [19]. We have verified that different choices of a, 
(3, 7 and 5 are consistent (see Appendix B). 

Once we have determined that the matching equation is integrable, we can proceed 
to solve it to find the coordinate transformation. Solving Eq. (34) we find 

1 



(X 



a 2 



1 



A a px? + 



2b 



3b 2 



[x k x k ] 



1 

If- 



1 rr> ry I ry /y . /y /y . 

, Ju JU fa J Ju Ju 2 <JO •* J % 



3 / 2 - 3x 2 



X k Xk 



tr 



^ " ^W)i{F Sl )iX d + (F af3 ) 2 xP + (C«) 2 



(40) 



as we show in detail in |Appendix B 



The quantities (F a p)i, (F a p) 2 and {C a )2 are determined when one carries out 
matching to higher order [19]. It should be noted, however, that since the IZ metric we 
employ here is in a different gauge than that used in [19] , the matching transformation 
is also different. Thus, one cannot simply use the results of [19] here. For simplicity, we 
will set (F a/3 ) 2 



a 2 



in the next section, while (F a p)i is determined in section 3.4 



3.3. Matching in the Spinning Case 

Let us now concentrate on matching in the spinning case and begin by estimating 
the order at which the spinning contributions would enter the matching calculations. 
Spin terms first enter the NZ metric at 0(v 5 ), 0{v A ) and 0(v 5 ) in the goo, g^ and g^ 
components, respectively. Spin terms first enter the IZ metric at 0(MA 2 /R 3 ) = 0(v 6 ), 
0{MA/R 2 ) and 0[M{VM 2 - A 2 - M)/R 2 ] = 0{v 4 ), and 0(MA 2 /R 3 ) = 0{v 6 ) in the 
9oo, Qoi and gij components, respectively. Here, since A has dimensions of mass, A/R 
is at most of 0(M/R) = Oiy 2 ). Also, the leading-order part of the coupling between 
z R/i, m and the BH spin arises at 0(AR ZR/i^ m ) = C(f 4 ), since ZR/i^ m oc m 2 /b 3 . 

This order counting argument suggests that spin contributions will first enter the 
matching calculation at 0[(m 2 /b) 2 } = 0(v 4 ). Therefore, one could take the spinning IZ 
metric and the spinning NZ metric, and apply the non-spinning matching transformation 
to obtain a spacetime that is properly asymptotically matched, without having to modify 
the matching transformations with spin terms. This is one the main results of this paper. 
We should note though that for many applications, the IZ metric must be written 
in appropriate harmonic, and horizon-penetrating coordinates, while the perturbation 
should be constructed with the appropriate boundary conditions [27] . 
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3.4- Restoring Temporal Dependence of the Tidal Fields 

During matching, all fields are expanded in the BZ, and thus, quantities that are time- 
dependent, such as the tidal fields, are Taylor expanded about t = m t/b 1. 
However, one can restore the time-dependence of these tidal fields, as discussed in 
Appendix B of [19] via the replacements 

Xi — > Xi cos ut + \ji sin cot , iji — )■ — Xi sin ut + y-i cos ut . (41) 

Using the above substitution, the zeroth order coordinate transformation becomes 

(X a ) = x a -b (x a cos ut + y a sin ut) 

m 



= x a -—bx a -J T ^J — ty a + 0(v' 2 ). (42) 
m y b \ m 

One can think of zeroth-order term, while y a is a first-order term. It is noted 



that the latter is of the same form as Eq. (30) if we rewrite it as - \J m 2 jmty a = 
\/m 2 /mip y a x^ where i = —1. We then have 

(F aP ) 1 = 2j—i v y a]i (43) 

where we used the antisymmetrization properties of the (F a p)i matrix that is necessary 
to leave the Minkowski metric in the zeroth-order matching unchanged in this coordinate 
transformation. As noted in [19], (F a/3 )i represents a boost (see also the low-order 
matching result of [20]). Therefore, we may derive the first-order time coordinate 
transformation from a Lorentz boost. 



Using the substitutions in Eq. (41), the quadrupolar field becomes 



2m 2 6m 2 6m 2 

ZRfl = -p- , z R<2 = -p- cos 2ut , z R -2 = -p- cos 2ut , 

6m 2 6m 2 
zi2 = — 7^-sin2a;t, zj _ 2 = — ^- sin2wt , (44) 
tr b 6 

in terms of ZRji )m and the other components vanish. 
4. Numerical Analysis 

Using the matched metrics of the previous sections, and the transition functions 
discussed in Appendix C , we can construct an approximate global metric. In this section, 



we study equal-mass BBHs and focus on the IZ and NZ only. The FZ metric becomes 
non-negligible at field points r/m > 40 and 110 for orbital separations of b = 10m and 



20m (see Eq. (C.5)), and it is automatically asymptotically matched to the NZ one by 
construction. 

Let us first look at the volume element of the 4-metric for a BBH with spins 
Xi = 0.9 = X2, where the dimensionless spin parameters xa = WaI/^a- Although the 
volume element is coordinate dependent, this is still a useful quantity to study when 
matching to verify that indeed the metrics approach each other smoothly in the BZ 
and in the given coordinate system. Figure [2] plots the volume element as a function 
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of x for different values of z and y = 0. This figure shows that indeed the IZ and NZ 
metrics smoothly match onto each other. The smooth matching exhibited by the volume 
element is characteristic of all metric components. 
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Figure 2. Global y/—g for an equal-mass BBH with equal and aligned spins 
(XI7X2) = (0.9,0.9). The orbital separation is b = 10m, and each BH is located 
on the x-axis at ±5m. The volume element is plotted as a function of x, with y = 
and different colors and styles correspond to different z-values. Since the figure is 
symmetric about x = 0, we only present results around BH1. The vertical dashed 
and dotted lines roughly correspond to the inner and outer boundaries of the IZ /NZ 
transition, respectively when z/m = 0. In the region outside the vertical dotted lines, 
the global metric reduces essentially to the NZ one, while in the region inside the 
dashed lines around x/m = 5 the global metric reduces to the IZ one. Elsewhere, the 
global metric transitions between IZ and NZ metrics. 

Let us now study how the matching behaves as a function of spin parameter. For 
this, we concentrate on the g$ y component of the global metric, as this is one of the 
components most affected by different spin values. The left panel of Figure [3] plots this 
component along the x-axis, for an equal-mass BBH with different spin configurations: 
(X1.X2) = (0,0), (xi,Xa) = (0.9,0), and (xi.Xz) = (0.9,0.9). We observe the strong 
spin-dependence of this metric component close to either BH (around x/m = ±5). 
In spite of this strong spin-dependence, we observe that the matched metric smoothly 
transitions between the IZs and the NZ in the BZs. Such a smooth transition is also 
shown on the right-panel of this figure, where we present both the global and the NZ 
only metrics for the (X11X2) — (0.9,0.9) case. 
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35 
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Figure 3. Left: Global go y along the ir-axis for an equal-mass BBH with different 
spins. The inset zooms in to the middle region of the transition from the IZ to the NZ. 
The IZ1 metric contributes at ~ 2% to the global metric x/m — 15. Right: Global and 
NZ only go y along the a;-axis for an equal-mass BBH with spins (Xi,X2) = (0.9,0.9). 
The vertical lines are the same as in Figure [2] 




Figure 4. Left: Ricci scalar R calculated with the global metric along the x-axis for 
an equal-mass BBH with orbital separation b — 10m but different spins. The insets 
zoom to a region close to BH1 for the (xi,X2) = (0,0) and (xi,X2) = (0.9,0.9) cases. 
In the insets, the dashed vertical lines denote the location of the event horizon of the IZ 
perturbed Kerr metric at x = (5 ± 0.5)m, and we stop plotting the Ricci scalar inside 
this region. Right: Same as left panel but for the (xi,X2) = (0.9,0.9) case plotted for 
different orbital separations: b = 10m and 20m. The x-axis is the field point distance, 
normalized to the orbital separation. 

Finally, let us use the Ricci scalar R as a measure of the accuracy to which the global 
metric satisfies the vacuum Einstein equations. An exact solution would of course satisfy 
R = 0, but since the metrics we are employing are approximate, their associated Ricci 
scalars will not vanish exactly. The questions one wish to address are the following: (i) 
are the IZ and NZ metrics as accurate as they are supposed to be? (ii) do the transition 
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functions introduce an error larger than that intrinsically contained in the approximate 
IZ and NZ metrics? Question (i) should clearly be answered in the affirmative, if the IZ 
and NZ metrics have been correctly calculated. Question (ii) should also be answered 
in the affirmative, provided the transition functions are built to satisfy the Frankenstein 
theorems |18j . 

The calculation of the Ricci scalar is not trivial due to the complexity of the IZ 
metric, the matching parameter and coordinate transformation, and the use of transition 
functions. For these reasons, we compute the Ricci on the t = spatial hypersurface 
through numerical methods. In doing so, we not only need the spatial derivative of the 
global metric, but also its time derivative. The latter requires knowledge of the time 
evolution of the orbital phase and frequency, summarized in the Appendix of [56], as 



well as the location of each BH in Eq. (15). We are using the currently-known highest 
PN results for the time evolution, i.e., the 3.5PN order for the non-spinning cases, and 
some lower order terms for the spinning cases [57]. When computing derivatives, it is 
of course critical to use sufficient numerical accuracy, as otherwise the Ricci scalar will 
be highly inaccurate, specially in the region close to the BHs. 

The left panel of Figure [4] presents the Ricci scalar for different spin values, while 
the right panel shows this scalar for different orbital separations. First, we see that 
the metrics are indeed approximate solutions to the Einstein equations, as evidenced 
by the right panel. That is, as the orbital separation is increased, the accuracy of the 
metrics also increases and the Ricci scalar decreases at the expected rate. Second, we 
see that the transition functions have been built to satisfy the Frankenstein theorems, 
as they clearly do not introduce errors larger than those inherently contained in the 
approximate metrics. That is, the height of the humps shown in the left panel of 
Figure [4] for x E [—15, —6] U [—4,4] U [6, 15] is of the same order as the error naturally 
inherent in the approximations. All of these conclusions are mildly sensitive to the spin 
values, demonstrating the validity of the global metric for all spins. 

The extreme Kerr limit, i.e., as xa 1, must be treated carefully. From the 
matching stand-point, this limit is perfectly well-behaved in a mathematical analysis 
sense. From a numerical standpoint, however, this limit is difficult because of excision. 
When evolving quantities numerically, one sometimes excises the region interior to the 
BH horizons, placing the excision boundary somewhere well inside the horizon. In the 
extreme Kerr case, the horizon shrinks and this could potentially lead to numerical 
problems with the excision boundary. We stress however that this is a numerical 
problem, and not a mathematical one with asymptotic matching. 



5. Discussion 



We constructed a global approximate metric that represents the binary inspiral of 
spinning compact objects. We split the spacetime into different zones: the IZs (close 
to either BH), the NZ (far from either BH but less than a GW wavelength from 
the center of mass) and the FZ (farther than a GW wavelength from the center 
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of mass). In each of these zones, the spacetime is approximated through either 
BH perturbation theory techniques in the IZs or PN techniques in the FZs. These 
approximate metrics are then related to each other via asymptotic matching, which 
provides a coordinate and parameter transformation that renders adjacent metrics 
asymptotic to each other in their mutually overlapping regions of validity. Once the 
metrics have been asymptotically matched, they can be stitched together via certain 
transition functions to yield an approximate metric for non-spinning BBHs [121 H7J HH] ■ 
Here, we have extended this work to spinning BBHs where the IZ metrics are modeled 
through perturbed Kerr spacetimes [27] and the NZ and FZ metrics include spin-orbit 
terms from PN theory [2J1 [221 IM] • After matching, the IZ metric reduces exactly to the 
Kerr solution in the limit of infinite binary separation. 

The approximate global metric constructed in this way is technically valid on an 
initial t = 0, spatial hypersurface, but it can be extended to capture the entire temporal 
evolution of the binary. This evolution will be described in a future paper [29J. Of course, 
this time-dependent, approximate global metric is only valid provided the approximate 
metrics in each zone remain properly asymptotically matched, which in turn holds true 
if and only if a BZ exists. These overlapping regions of validity are dynamically squeezed 
out as the binary shrinks, and we expect the global construction to fail for sufficiently 
small separations, as will be studied in [2H]- When the BHs are close enough, numerical 
relativity is required to describe the BBH spacetime. 

We are planning to extend this work in multiple ways in the near future. We have 
here carried out asymptotic matching to lowest order, and thus, a clear extension would 
be to re-do this calculation but to higher order. For this to be feasible, however, one 
would first have to construct a vacuum-perturbed Kerr metric that includes the first 
time derivative of the electric and magnetic quadrupole tensors, as well as the octupole 
tensors. Such a metric is not available at this time, which is why we were forced to stop 
the matching at lowest order. 

Another avenue of future work would be to explore how the dynamical approximate 
global metric constructed here affects certain astrophysical phenomena, such as accretion 
disks around a BBH. A similar study was carried out in [I], except that there non- 
spinning BHs were considered in the NZ metric. We expect spin-orbit coupling to be 
important to properly account for the total angular momentum budget of the system. 
For example, the spin-orbit coupling can have a dramatic effect on the inspiral rate. 
When the BH spins are aligned (or even partially aligned) with the orbital angular 
momentum, the merger is delayed, while when they are anti-aligned, the merger happens 
much more quickly. This effect, also know as "hangup" effect [SB], is also responsible for 
very large kicks of the final merger remnant [59]. We also expect spins to be important 
to properly describe the dynamics of the individual accretion disks that may build-up 
around each BH during the inspiral process, specially at large separations. It would be 
interesting to examine what the differences are when spin is included, and whether this 
would lead to an electromagnetic observable precursor to BBH mergers. We will explore 
this in a forthcoming paper. 
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Appendix A. Weyl Scalar for the IZ and NZ Metric 

Our starting point is the IZ metric, which is given by Eqs. (YG-45) and (YG-46), where 
Eq. (YG-n) denotes Eq. (n) in [27]. Since the m = mode is of higher-order in the 
analysis of [27J, it was ignored there. In this paper, however, we wish to include this 
mode, and thus, we derive it following the methods in [27] . 

First, the radial function R m in Eq. (YG-8) is rather simple when m = 0: 

R = const. = 1 , (A.l) 

where this constant is determined by requiring that R m — > 1 for large radius. One can 
also think of this solution as that for the m/fl modes in the Schwarzschild limit |41j . 
Using Eq. (YG-21) for R , we calculate the potential The radial dependence is simply 
A 2 , and we can directly extend Eq. (YG-34) to including the m = mode. Thus, 
although there is a singular coefficient in the Schwarzschild limit (see Eq. (YG-33)), 
we may use Eq. (YG-34) for the non-spinning casdjj] If we considered time-dependent 
perturbations, the reconstruction of the metric would be more complicated, as discussed 
for example, in [391 HP] , 

Next, using Eq. (YG-45) with all (£,m) = (2,m) modes, we calculate the Weyl 
scalar ipQ. One finds that ipo = — (l/2)ip^ ITlg ^ where ^Q° ng ' ) denotes the Weyl scalar in 
Eq. (YG-5). Therefore, we have obtained a relation between £ kh B^i and z R /^ m , as given 
in Eq. (22). The factor of 2 difference arises due to a difference in normalization of ipo 



and \1/ (see, e.g., (15) in [61]). In practice, the time-time and time-space components 



shown in [27J and [T5] are the same in the M/R — > limit with Eq. (22). 

The Weyl scalar ipo can be calculated directly from the NZ metric, independently 
from section [3] Calculating the Weyl tensor from the NZ metric and deriving the Weyl 
scalar by contracting with the tetrad, we can compare it with the IZ ipo in the BZ. From 



Eq. (21) for the NZ metric in the BZ, we have 



^ = [(1 + cos 2 9) cos 2(p - sin 2 9 - 2 % cos 6 sin 2<p] , (A.2) 
2 o 6 

|| In the slow-motion approximation, i.e., when the characteristic velocity »<1, the time dependence 
of zi m in the potential "J is also slow, and we do not need to consider a rapidly spinning BH [27| . 
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to leading order in 777,2/6, where we are using spherical polar coordinates and we ignored 
any coordinate difference between the IZ and NZ. Matching then forces z^q = 2rri2/b 3 , 
zr,±2 = 67772 /b 3 and all other coefficients to vanish. The above result is converted to 



£ij = (m 2 /b 3 )(5ij - 3xiXj) via Eq. (22). 



Appendix B. Detailed Verification of the Matching Calculation 

We have so far studied matching for a given set of indices in the integrability condition 



of Eq. (35), but now we extend this result to all indices. Using the same notations as 



before, and taking (a, j3, 7) = k) 7^ and 5 = 0, we now have 

^•(^0)2-^(^0)2 = 0, (B-l) 
which is explicitly written as 

= — Sjk , D t 3 (4m)o(^ n )o(^ m )o + 25jfc-— r (£in)o(X n ) — 2 ■jj— I (£ in ) (X n ) 
[ K )o \ K )o Wo 

— $ij / {£nm)o(X U ) (X m ) + 3 -. J * (£ nm ) (X n ) (X m ) 

— $ik , D { 3 (4m)o(^")o(^ m )o + 25 ife (£j n ) (X n ) — 2 % \ (Sj n )o(X n ) 
+ 2 7^~(^i)o + {£ k n)o{X n )o + 2—^—(£ ki ) + 2-^—(£ kj ) 



-27M(Uo(n-(^^ (B.2) 
\ n )o 

where x % = (X l ) . This equation is trivially consistent, and thus, it provides no 
additional information about (£ij)o- 

The last set of indices to verify is (a, /3, 7, 5) = (i,j,k,l) 7^ 0. The integrability 
condition becomes 

dijiSuh + duiSyh - d u (S kj ) 2 - d kj (S u ) 2 = . (B.3) 

As before, this equation can be divided into a polynomial and a non-polynomial part. 

The non polynomial part gives again (iWi) = mi, and once more the constant and first 

order pieces give trivial relations. The quadrupolar part gives 

2 - 4 2 - 

= (6xiXj - 26 i:j )8 k i + -(£ij) S k i + + (Qx k xi - 26^)5^ + -(£ k i) Sij 

2 

- (6x k Xj - 28 kj )8u + -(£ kj )aSu , (B.4) 



which once again leads to (£y)o = <5ij — 3x 



Let us now focus on the coordinate transformation. We need to solve 

2(A M) ) 2 = G$) 2 - (^) 2 - (F a 7 ) 1(^7)1 , (B.5) 
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or more explicitly, using the expression for (£kj)o, 



2(A M ))2 = [(2 - ^x)A a p - (Fa')i(F^)i 



, 2 / _ 2 



+ 



1 / / 00 ■ 



(B.6) 



where we used the abbreviation 3/ OC 00 00 OC £y because xo = 0. The solution to this 
equation is obtained by adding the general flat solution to a particular solution. The 
general solution is then 

(X a ) 2tg = {F aP ) 2 x p + {C a ) 2 . (B.7) 



The first bracket in Eq. (B.6) is the same as the one found in [T9] during second- 
order matching. The particular solution for this term is 

A^x 5 . 1 

X„. 1 

2 



(X 



1 



x 



2b 



For the second bracket, a particular solution is 

For the third bracket, we can take the particular solution 

1 



(B.8) 



' ry rf> \ rjr> rjr> , T* T* ■ 

i JU JU I »A/ »A/ j «/y JU -fa 



SI. 



{X, 



a)2,p3 



3b 2 



oc oc 3 oc , \f oo oc 



x k 



(B.9) 



(B.10) 



Combining the general and particular solutions, we obtain the expression in Eq. (40). 



Appendix C. Transition function 

The construction of a smooth approximate global metric requires the stitching of the 
asymptotically matched IZ, NZ and FZ metrics through certain transition functions, 

fla.ii /near; /inner, 1 and /inner, 2- 

9fJ.u = (1 /far) "I /near [/inner, 1 9^ u ~H (1 /inner, l) 9)iv ^] 

+ (1 — /near) [/inner,2 9^1 + (1 ~ /inner,2) fi^f^] j + /far 9^ 

These transition functions can be modeled via 
' 0, r < r , 



7 FZ 



(ci) 



/(r, r , w, g, s) 



1 + tanh 



- (x(r, r , w) — r 

7T V X{r, r , w) 



(C.2) 



r < r < r + w, 
1, r > r + w, 

where x( r , r o, w ) — tan[7r(r — r )/(2w)], and r , w, q and s are parameters. Such 
a transition function has been greatly discussed in [15], [T71 [19] and it satisfies the 
Frankenstein conditions of |18|. 
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The transition functions are chosen with the following parameters. 

/far = f(r, A/5, A, 1, 1.4) , (C.3) 

/ near =f(x, 2.2m 2 - — , 6 -2.2m, 1, 2.5), (C.4) 

m 

/i„ner,A = f(r A , 0.4r5, 3.5 (m& 4 ) 1/5 , 0.2, 6/m) . (C.5) 

Here, r is the distance from the binary's center-of-mass to the field point. The parameter 
s = 1.4 has been discussed in [29] and it is different from that of [19J, which reduces the 
magnitude of the first radial derivative. The quantity A = ix^fWjM is the gravitational 
wavelength in the Newtonian limit. We assume that the BHs are initially located on 
the x-axis. 

The transition function fi nnert A is somewhat different from that chosen in [TJ] . This 
is because we match here to a different order than in the latter paper. In particular, 
the "transition radius" r T Al where the leading-order uncontrolled remainder of adjacent 
zones becomes comparable, is here given by 




(C.6) 



where the left- and right-hand sides correspond to the IZ and NZ remainders, 
respectively. The quantity {r^/b) 3 in the IZ remainder arises due to ignorance in the 
octupole contribution to the BH perturbation. In the NZ, the remainder is (m^/r^) 2 . 
The transition radius is thus r T A = (m^6 4 /m) 1 ' 5 , which is the same as in [T7J. The 
numerical coefficients in /i nne r,A, 0.4 and 3.5, are also from [17] . and the (mfr 4 ) 1 / 5 
dependence is obtained by setting = m in r A , as studied in [T9~] . 



Appendix D. Higher Order Metrics 



Ref. [19] constructed a 2.5PN, non-spinning NZ metric following [23]. In this paper, 
we have included certain higher PN order terms, when numerically calculating certain 
quantities in section [4] 

For the NZ metric, this higher-order extension is achieved by adding the following 
terms, 5g™l to the 1.5PN order metric of Eq. (pi): 
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In the above equations, we have introduced the notation 
S = ri + r 2 + b , 



(D.2) 
(D.3) 



where the quantity S is a distance parameter that is not to be confused with the 
magnitude of the spin angular momentum s% or s 2 . The other higher-order pieces can 
be obtained up from Eq. (7.2) in [23]; for example, Sg^Q is a contribution of 0(1/ c 6 ) 
that can be found in Eq. (7.2a) of [23J. 



with 



For the FZ metric, the higher-order extension is obtained by replacing Eq. (12) 
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with 0(i; 6 ) remainders. The metric potentials hp^ must also extended to higher-order 
via 
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again with 0(v 6 ) remainders. 
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